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In this paper, we report the applicability of the density matrix renormalization group (DMRG) 
approach to the cylindrical single wall carbon nanotube (SWCN) for purpose of its correlation 
effect. By applying the DMRG approach to the t+U+V model, with t and V being the hopping 
and Coulomb energies between the nearest neighboring sites, respectively, and U the onsite Coulomb 
energy, we calculate the phase diagram for the SWCN with chiral numbers (m = 3,ri2 = 2), which 
reflects the competition between the correlation energy U and V. Within reasonable parameter 
ranges, we investigate possible correlated groundstates, the lowest excitations and the corresponding 
correlation functions in which the connection with the excitonic insulator is particularly addressed. 
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Since the discovery of carbon nanotube, many ef- 
forts, both experimental and theoretical, have been 
put into understanding the electron correlation ef- 
fect, specially for thin single wall carbon nan- 
otubes(SWCN). The transport as well as the angle- 
integrated photoemissioni*2studies on the bundles of 
SWCN show interesting evidence of the non-trivial cor- 
relation effects which are believed to be contributed 
by the metallic SWCNs and are consistent with the 
Tomonaga-Luttinger theorji 3 ^^. For the semiconduct- 
ing SWCNs, the discrepancy between the experimen- 
tally measured ratio of the absorption frequency to 
the emission frequenc}*^ and that predicted by the 
tight binding approximation(TBA) 8 is also attributed 
to the correlation effects in Iiteratures2*i2*ii. Theoret- 
ically, the correlation effect of SWCNs so far has been 
treated^^iiSiiiiiSiiiii by ab initio numerical calcu- 
lations, perturbative renormalization group(RG), and 
various improved or modified Hartree-Fock approxima- 
tions (HFA). However in these approaches the strong cor- 
relation effects are not easy to be treated seriously. 

It is known that the 7r-electron has its wave function 
extending perpendicular to the surface of SWCN 8 with 
a weak overlap among themselves, so that its kinetic en- 
ergy in a range of 2.4-3.2eV is much smaller than the 
onsite Coulomb energy about 11.76e\4i^ while the near- 
est neighboring(NN) Coulomb energy is also as big as 
an order of 5eV. Therefore the SWCN is a system which 
might exhibit strong correlation effect. We consider an 
extended Hubbard model to study the correlation physics 
of SWCNs 

H = 2_j C i,<r C J,a + h - c - + U ^2 ni <T n U 
<i,j> i 

+ vj^im-i)^-!), (l) 

where the hopping energy t is set to be unit, while U and 
V are onsite and NN Coulomb energies, respectively. Cj 



and c] are the electron annihilation and creation oper- 
ators with i the site index, respectively. The particle 
number operator for spin a at site i is n,- lCT = cj a Ci t(J 
and rii = + n^. Here we neglect the curvature 
induced difference among the three hopping directions 
on each site. This minimal model has been studied in 
Refs.|3IT^| for metallic tubes by perturbative RG and un- 
restricted HFA, respectively. It has also been applied to 
the graphite system^. 

In this paper, we report the applicability of the DMRG 
approach to the cylindrical SWCN for purpose of study- 
ing its correlation effect. By applying this approach to 
the t+U+V model Eq.(l), we compute the phase dia- 
gram for the SWCN with chiral numbers (ni = 3,ri2 — 
2), which essentially reflects the competition between the 
correlation energies U and V. Within reasonable pa- 
rameter ranges of U/t and V/t, we investigate possible 
correlated groundstates, the lowest excitations and the 
corresponding correlation functions. 

As well accepted, the DMRG approach is an accu- 
rate and efficient treatment for the quasi-lD strongly 
correlated systems^. The main problem for the imple- 
mentation of DMRG to the SWCN lies in how to con- 
struct an appropriate superblock configuration for such 
a cylindrically warped hexagonal lattice sheet. We no- 
tice that, following the literature [lg, for a given pair of 
chiral numbers (n-i, ri2) with N the greatest common di- 
visor, the corresponding SWCN can be constructed by 
successive piling of motifs along the tube axis conducted 
by a screw operation 18 :, where each motif contains N unit 
cells in coincidence with the TV-fold rotational symmetry 
of the tube. This implies an iV-multiple helical struc- 
ture hidden in the SWCN. The unit cells are labelled 
by two integers ijn,l) with I = 1,2, ...,7V which corre- 
sponds the l-th cell in the m-th motif. The two car- 
bon atoms in one unit cell are distinguished by an index 
v = A, B. Based upon the above observation, we map 
such TV-multiple helical SWCN onto a 2N chain lattice 
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so as to make the DMRG procedure can be straight- 
forwardly implemented. In this helical coordinate sys- 
tem, the site (m, I, A) has NNs labelled by (m,l,B), 
(m-^,l + pi,B), (to + ^,l-p 2 ,B), where pi,p 2 are 
integer solutions of equation n\p 2 — n 2 p\ = In par- 
ticular, for case of N = 1, which is the most advantaged, 
all the unit cells on the tube can be threaded by only 
one helical curve and I becomes a dummy index. Then 
the Hamiltonian in Eq.(l) can be rigorously mapped to a 
twisted two-chain model with finite long range hoppings 
as H = B t + H v + By with 

Ht = C m^ C m+^ + k.C. 

B V = U Y n m,f n m,l 
m,v=A,B 

By = V Y (»m " " !) ( 2 ) 

where S = {0, -n l7 n 2 }, n v m a = c v ^ a c v m al n v m = n v mA + 
n m I anc ^ c m,cr i s exactly the same operator ci. a but with 
its index relabelled in terms of motif-cell labels. Fig. la 
and lb show such a mapping for the tube (n\ = 3, n 2 = 2) 
with N = 1. 




FIG. 1: Illustration for the mapping from a cylindrical 
SWCN(la) with index (3,2) to a twisted two chain model 
with finite long range hopping (lb), where the solid lines con- 
nect the NN pairs. For purpose of DMRG calculation, we 
further squeeze the two chains in (lb) into a single chain (lc) 
with A atom and B atom arranged alternatively, where the 
lines are kept with the same meaning as in (lb). 

To extract most physical implications with small 
enough computing effort, in this paper, we apply the 
standard DMRG approach to the intrinsic semiconduct- 
ing SWCN (3, 2) with N = 1, in which the open boundary 
condition(OBC) is engaged. The computed groundstate 
phase diagram is plotted in Fig. 2. It can be divided into 
three phases as the excitonic insulator (E), normal semi- 
conductor^), and Mott insulator(M), as will be clear 
shortly. 

For our minimal model as Eq.(2) when both U and 
V are small, the system is a semiconductor which can 



be divided into two regions N s and Nt according to the 
first excitations. The boundary between them is shown 
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FIG. 2: The phase diagram of semiconducting SWCN(3,2). 
See the text for details. 

as the dashed line in Fig. 2, across which a level crossing 
takes place for the gapped excited states with different 
spin symmetries, i.e., the singlet in the N s region and 
the triplet in the Nt region, which are favored by V and 
U, respectively. We interpret these excitations in con- 
nection with the singlet and triplet excitons for the in- 
trinsic semiconducting SWCN. The groundstate in these 
two regions is a non-degenerate singlet state, of which the 
equal time correlation function(ETCF) is exponentially 
decaying in both spin and charge sectors. The difference 
is that density-density correlation is much stronger than 
the spin-spin correlation in the N s region and vice versa 
in the Nt region. 

If we still keep U small but increase V, the first spin 
singlet excitation would get lowered. When it reaches the 
groundstate, the system then undergoes a transition into 
a new phase in association with the excitonic insulator 19 . 
As a consequence of the "fallen state" , there are two mu- 
tually orthogonal degenerate groundstates in the E phase 
as \G + ) and \G_ ), in which the 2-fold degeneracy is 
understood in the thermodynamic limit. Both of them 
are spin singlet and eigenstates of the parity, belonging to 
eigenvalues ±1, respectively. Moreover, the amplitude of 
the density-density correlation increases with increasing 
V as shown in the upper panel of Fig. 3, until into the E 
phase it is identical for both the two degenerate ground- 
states and no longer decays. This non-decaying corre- 
lation could be interpreted as the existence of a charge 
density wave in the two groundstates. When comput- 
ing the matrix element of the particle number operator 
n v m between these two degenerate states, it is found that 

the diagonal term (G^Vml^) {or (G^ln^G^)) is 
trivially equal to 1, but the off-diagonal matrix element 
is 

(G? ) K|G^ ) )~(-l)"n (3) 

with n being a site-independent constant and (—1)" = 
±1 for v = A, B, respectively. By mixing the two states 
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one can construct the groundstate with alterna- 
tive density distribution in the A and B sublattices. This 
can be viewed as the exciton(particle-hole pair) conden- 
sation in the groundstate of E phase with a spontaneous 
symmetry breaking between the sublattices. 



excitation gap is opened. It decreases with increasing 
U and vanishes at about U — 8.8. Fig. 4b with V — 
0.8 shows a different picture. The singlet gap emerges 
gradually from zero at the E — N s phase boundary. When 



^o— V=0.7 — a— V=0.8 — V=0.9 — o— V=1.0| 

S222S222S2222$$2222??ooooo<> 





© U=1.0- 



♦Iff* 



V=2.0. 



— •— U=5.0 — U=6.0 
U=8.0 — U=20.0 



10 



20 



30 



40 



50 



60 




FIG. 4: The first excitation gap as a function of onsite 
Coulomb potential U, for V = 2.0 in (4a), V = 0.8 in (4b), 
and V = 0.5 in (4c). The vertical dashed lines indicate the 
transition points. 



FIG. 3: The charge and spin correlation functions. The coor- 
dinate system for the horizontal axis is explained in Fig. lc. 
The origin is set on the atom A in the middle cell of the tube. 
The data are obtained by keeping 2000 states and with 1-2 
finite size sweep to get convergent results. The upper panel 
gives the change of charge correlation C(|i — j\) with respect 
to V for fixed U = 1.0 and the lower panel is the antiferro- 
magnetic spin correlation S(\i — j\) for different U with fixed 
V = 2.0. 

Starting from the region Nt, if we further increase U, 
the triplet first excitation would be lowered continuously 
with the groundstate spin-spin ETCF growing. In the 
lower panel of Fig. 3, we show the variation of spin cor- 
relation versus U for fixed value of V. Until some critical 
value of U = U c , the triplet spin gap vanishes in the 
thermodynamic limftsSi, and a transition from phase N 
to the Mott phase M takes place. In the Mott phase, 
besides a finite charge gap the system develops gapless 
spin wave type excitations. The Mott phase is far from 
the possible realistic parameter region of U and V for the 
SWCN(3,2), however it may be realized in the small gap 
semiconductor or metallic tube. 

The phase boundaries in Fig. 2 are scanned by comput- 
ing the energy gap in the particle-hole channel for various 
U and V. In this procedure, we adopt the infinite algo- 
rithm with 1000 states being kept and then extrapolate 
the finite size data to the thermodynamic limit. The cal- 
culated phase boundaries will also be checked by other 
independent method below. In Fig. 4, we give the vari- 
ation of the gap with respect to U for some fixed values 
of V as examples. 

For V — 2 in Fig. 4a, the 2-fold degeneracy of the 
groundstates in E phase is lifted abruptly at the E — 
Nt phase boundary Uc = 5.6 and a finite spin triplet 



the increasing U reaches the N s — Nt boundary, a level 
crossing takes place from the singlet excitation to the 
triplet excitation. Fig. 4c with V = 0.5 only shows a iV s — 
Nt transition qualitatively the same as that in Fig. 4b. 

To further explore the nature of the above shown phase 
boundaries, and in particular, as an independent check, 
we compute the average of double occupancy and the 
NN density correlation over the groundstate, denoted by 
6 and p, respectively, as in Ref.|20| 



6 = dE GS /dU = (Hu/UL) 
p = dE GS /dV = (Hy/VL) 



(4) 



Here Eqs is the groundstate energy per site and the 
Hellmann-Feynmann theorem is applied. For fixed V = 
2.0 or 0.8, the 8 and p as functions of U are plotted in 
Fig. 5. The calculations are performed on the tube with 
120, 160 and 200 carbon atoms by keeping 2000 states 
and sweeping 1-2 times to get convergent results. The 
results calculated from different atom numbers coincide 
with one another nicely, which provides a strong evidence 
that the calculated main features of the transition will 
survive in the thermodynamic limit. 

As shown in the right panel in Fig. 5, there is an abrupt 
jump for both 6 and p at the critical point V = 2.0 and 
U = 5.6. It indicates the phase boundary between E 
and Nt belongs to the first order transition. Note that 
the phase boundaries (dashed lines) in Fig. 5 are extracted 
from the phase diagram Fig. 2, which agree with the data 
of S and p very well. If we further increase U from region 
Nt to phase M, 5 and p vary smoothly across the phase 
boundary (see also right panel of Fig. 5). It is understood 
that the Nt — M phase transition is of the Kosterlitz- 
Thouless type without symmetry breaking, where the Nt 
region with a gapful (massive) triplet excitation spectrum 
transits into the M phase with gapless (massless) spin ex- 
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citations. In the left panel of Fig. 5, there is a disconti- 
nuity in the slope of both S and p as functions of U on 
the N s — E phase boundary, which means this bound- 
ary belongs to a second order phase transition. Across 
the N s — Nt division line, a level crossing of the excited 
states, not a phase transition, takes place. 
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FIG. 5: The plot of 8 and p as functions of U for fixed V. The 
vertical dashed lines indicate the transition points which are 
extracted from the phase diagram Fig. 2. The calculated data 
are for 120 atoms(open squares), 160 atoms(open triangles), 
and 200 atoms(cross symbol), respectively. 



Briefly the transitions from the semiconductor into the 
excitonic insulator and Mott insulator take place when- 
ever the band gap is overwhelmed completely by by V 
and U, respectively. For other semiconducting SWCNs, 
the augment of the tube diameter will reduce the TBA 
band gap, which will consequently shift the phase bound- 
aries E — N s and Nt — M. However the phase boundaries 
E — Nt and N s — Nt should not change too much be- 
cause they can be viewed as the balance position of the 
competition between U and V. 

In our minimal model, the longer range(LR) part of 
Coulomb interaction is ignored, which is nonessential in 
describing the correlated groundstate of semiconducting 
SWCNs. Since we are dealing with the low energy exci- 
tations of the intrinsic semiconducting SWCN, the main 
effects of the LR Coulomb interaction are confined in 
the forming particle-hole excitations and its remnant ef- 
fect leads to a renormalization of U and V. Based on 
the above argument, the groundstate of SWCN may be 
an excitonic insulator, normal semiconductor with differ- 
ent type of first excitations and Mott insulator, which 
is dependent on the chiral structure of SWCN and also 
the environment. In experiments, people can realize the 
transition of the correlated groundstate of SWCNs by ap- 
plying a uniform pressure or by changing the band gap 
via the magnetic fhui 22 i 23 i 24 . 
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